DISTRIBUTIONS WITH DYNAMIC TEST FUNCTIONS AND 
MULTIPLICATION BY DISCONTINUOUS FUNCTIONS 
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Abstract. As follows from the Schwartz Impossibility Theorem, multiplication of two distri- 
butions is in general impossible. Nevertheless, often one needs to multiply a distribution by a 
discontinuous function, not by an arbitrary distribution. In the present paper we construct a 
space of distributions where the general operation of multiplication by a discontinuous function 
is defined, continuous, commutative, associative and for which the Leibniz product rule holds. 
In the new space of distributions, the classical delta-function 5^ extends to a family of delta- 
functions 5", dependent on the shape a. We show that the various known definitions of the 
product of the Heaviside function and the delta-function in the classical space of distributions 
V become particular cases of the multiplication in the new space of distributions, and provide 
the applications of the new space of distributions to the ordinary differential equations which 
arise in optimal control theory. Also, we compare our approach of the Schwartz distribution 
theory with the approach of the Colombeau generalized functions algebra, where the general 
operation of multiplication of two distributions is defined. 



1. Introduction 

The theory of distributions created by L. Schwartz |Sch50| in the 1950s is an important part 
of modern mathematics. However, in 1962 R. Courant pointed out certain insufficiencies of the 
Schwartz distributions as generahzations of functions, in particular, the absence of a general mul- 
tiplication of distributions CH62| . This multiplication is important for the definition of a solu- 
tion for a wide class of ordinary and partial differential equations Cq185j . It was mentioned by 
J. F. Colombeau in |Col85| that there is a strong belief by many mathematicians that a general 
multiplication of distributions is impossible in the Schwartz distribution theory, unless we drop 
essential properties (such as continuity), which is not acceptable for applications to differential 
equations. 

Along with that, in many cases in the theory of differential equations one needs to be able to 
multiply a distribution not by an arbitrary distribution, but by a discontinuous function only (which 
is also a distribution). For instance, the product of a distribution and a discontinuous function 
arises in ordinary differential equations and partial differential equations of optimal control theory, 
e.g. see |MR96. SV9 7, .M;!!).'? Dcr88 , where a distribution represents an impulsive control. In our 
paper we construct a new space of distributions which are extensions of the classical distributions 
from the space of continuous test functions to the space of dynamic test functions. In this space the 
operation of multiplication of a distribution by a discontinuous function is defined, commutative, 
associative and continuous. 

Further we provide some applications of a new space of distributions to the theory of ordinary 
differential equations. Also, let us mention that since the operation of multiplication of a dis- 
tribution and a discontinuous function can be induced by the operation of multiplication of two 
distributions, its properties allow us to compare, in a certain sense, various definitions of the prod- 
uct of two distributions in the classical space of distributions with continuous test functions V , 
which were proposed, in particular, in |Bag95[ [Bag02| ISar95l ISar94l ISar(l3) . in attempt to obtain 
a better definition. 

The operation of multiplication of a distribution by a discontinuous function (as a particular case 
of a general operation of multiplication of two distributions) is defined in the Colombeau generalized 
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functions algebra, which was apphed successfully to many problems related to ordinary and partial 
differential equations jCol85, ,CLNP89, CH92, Col92, CH093, CH094, CM94]. In the present 
paper we show that the commutative, associative and continuous operation of multiplication of a 
distribution by a discontinuous function can be also defined in the classical setting of the Schwartz 
distribution theory. Below we compare our approach based on the Schwartz distribution theory 
with the approach of the Colombeau generalized functions algebra (Section 4). 

1.1. Continuity properties of the operation of multiplication. Let us describe some natural 
requirements imposed on continuity properties of the operation of multiplication. Below the term 
"regular distribution" refers to a distribution which corresponds to a locally-summable function. 

The operation of multiplication of a distribution by a discontinuous function has to have the 
following property: 

Property A. If g is a discontinuous function, v is a distribution, {wn}^i is a sequence of reg- 
ular distributions, Vn — > v, then gVn — > gv. If v is a regular distribution, then gv is regular and 
corresponds to the usual product of g and v. 

A space of distributions where the operation of multiplication by a discontinuous function has 
property A was constructed in 'DK05 . This operation of multiplication arises from the following 
ordinary differential equation [Mil93, . 

(1) f{t,x)+g{t,u{t))v 

where w is a distribution (impulsive control), /, g are continuously differentiable, u is a control, 
which is in general discontinuous function |Mil98| . 

However, property A is too weak for consideration in the following ordinary differential equation 
with distributions which arises in many problems of optimal control theory (SZ97. ■SV97.. ■MR96) . 

(2) ± = f{t,x)+g{t,x)v, 

where d is a distribution, /, g are continuously differentiable. In the general case, if the distribution 
V is not regular, it is natural to expect that a solution x is discontinuous |HZ97itSV97nMR96 Mil93]. 
Thus, the product of a distribution and a discontinuous function in (0) must have the following 
property. 

Property B. If g is a discontinuous function, v a distribution, {vn}'^^i a sequence of regular dis- 
tributions, Vn V, and {ffnl^^ sequence of continuous functions, g„ — > g, then gnVn — > gv. 
If V is a regular distribution, then gv is regular and corresponds to the usual product of g and v. 

Thus, in order to define the continuous operation of multiplication of distributions by discontinu- 
ous functions possessing property B, which is important for the applications to ordinary differential 
equations, one needs to distinguish different approximations {wn}^i of the distribution v in T>' , 
and different approximations {gn}^=i of the discontinuous function g. 

This idea was formalized in CLN P89I ICH92L ICol92| (see the references therein) in the frame- 
work of the Colombeau generalized fimctions theory (see more detailed discussion in Section 4) 
and, in less generality, in (Mil89k iMil9 3 BR94, SZ97j (see further references therein), where the 
consideration of the sequence {vn}^=i as a part of the system Q in order to provide the uniqueness 
of solution was proposed. 

In the present paper we formalize this idea in the classical setting of the Schwartz distribu- 
tion theory, i.e., in the space of distributions with dynamic test functions, where the operation of 
multiplication by discontinuous function has property B. 

1.2. Distributions with dynamic test functions. Let us briefly describe our approach. In 
distribution theory (as a theory of space T>') the product of a distribution and a continuous function 
is defined by 

(3) {gf,v)^{f,gv) 

where / is a distribution, 5 is a continuous function and ip is a continuous test function. If g is 
discontinuous, then g(p is no longer a test function, thus the value of the right-hand side of Q is 
undefined. In order to employ the definition we need to extend the space of test functions so 
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that it contains discontinuous functions. In this case the value of the right-hand side of (|3J| would 
be defined for discontinuous g also. 

In the present paper we construct the space of distributions with the test functions possibly dis- 
continuous. We show that any classical distribution can be extended from the space of continuous 
test functions to our space of test functions. In particular, the classical delta-function 5r extends 

1/2 

into a family of delta-functions {5"}, where the parameter a such that J_i/2 (^{t)dt — 1, is called 
the shape of a delta- function (5". In the new space of distributions the product of the Heaviside 
function 9^. and the delta-function 6" is given by 

(4) 0.5^ = l^j^ ' ait)dty^, 

1/2 

where a{t)dt e R is a constant, and f3 is another shape. Similarly, the family of delta-sequences 
extends into subfamilies of delta-sequences having different shapes. Then a result similar to l@J 
can be obtained if delta-functions in are replaced by terms of delta-sequences having the same 
shapes. 

We should mention that generalizes known definitions of the product of the Heaviside func- 
tion and the classical delta-function in V given by 

(5) OrS-r — c5r, 

where c G R, which were proposed (in particular, through the multiplication of two distributions) 
in |KB98l lFil88 Tvr02, SZ97I ISar95l ISar94l ISarllS] . in order to obtain the optimal (from a certain 
point of view) value of c S R, e.g. definition Q for c = 1/2 satisfies the formal Leibniz rule for 
the differentiation of the product. 

In contrast to definitions of form (jsj, the operation of multiplication in the new space of dis- 
tributions is continuous, commutative and associative, the operation of multiplication in the new 
space of distributions satisfies the Leibniz product rule. 

Apart of the problem of multiplication of distributions, our study is motivated by the problem of 
the correct definition of the solution of an ordinary differential equation with distributions. Further 
we provide an application of our space of distributions to nonlinear ordinary differential equations 
containing one delta- function S", which arise in optimal control theory jSZ97[ ICTR96I IS VgTI IMil93| . 

Note that the set of points of concentration of delta-functions in v and even its cardinality, as 
well as the set of points of discontinuity of x, can not be fixed a priori, since in many optimal 
control problems it is a subject of optimization |SV97I lMR.961 ISZ97j . This leads to the necessity of 
considering the space of functions of bounded variation, and, as a result, to certain complications 
in contrast to the case when the set of points of discontinuity of a function, which is multiplied by 
a distribution, is finite and fixed. 

In order to construct this space of distributions we introduce the supplementary notion of 
the dynamic function. Most of the definitions for the ordinary functions (such as limit at point, 
support, boundedness, continuity, bounded variation) are transferred to dynamic functions without 
significant changes. 



2. Dynamic test functions 



2.1. Notations. Let / = (a,5) C R be a fixed open interval, in general unbounded. By R" and 
M„ we denote the space of vectors and the space of square matrices of order n with real elements, 
respectively. We denote by T{g) C / the set of points of discontinuity of a function g : I ^ R. 
Consider the algebra of functions g : I ^ H possessing one-sided limits 

g{a+) = lim g{t), g{b-) = lim g{t), g{T+) = lim g{t), g{T-) = lini g{t) 

t^a+ t^b— t—>T+ t-^T — 

for any r G /. We identify functions having the same one-sided limits, and possibly having different 
values in their points of discontinuity. We denote this algebra of functions (i.e., equivalence classes) 
by G = G(/) and endow it with the norm 

llfflk = supmax{|g(i+)|, \g{t-)\}. 
tei 
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With this norm G is a Banach algebra ! Uer02| . The algebra G is called the algebra of regulated 
functions Die69]. The set of points of discontinuity T{g) = {t € I : g{t+) ^ g{t—)} of a regulated 
function g G G is at most countable |Die69| . 

For a function .g G G let us consider a partition d = C /, ti < • • • < t„ of the interval /. 

Denote 



(total variation) . If vaTi{g) < oo, we say that 5 is a function of bounded variation on /. We denote 
this algebra by BV = BV(/), and define the norm 

llffllBV = \g{a+)\ +va.ii{g). 

It follows that BV is a Banach algebra. 

We denote by C C(/) the subalgebra of continuous elements of G. Denote by CBV the 
subalgebra of continuous elements of BV. Denote by L = ^{I) the Banach algebra of functions 
Lebesgue summable on / with the integral norm. Further AC = AC(/) stands for the Banach 
algebra of absolutely continuous functions DS62 . Let L(loc) and AC(loc) be the algebras of 
locally summable and locally absolutely continuous functions, respectively. 

We denote crt{g) = g{t+) — g{t—)- 

Let us denote by F = F(/) the algebra of functions / ^ R with the operations of point-wise 
addition, multiplication by an clement of R and multiplication. The elements of F are called the 
ordinary functions. 

The definitions given above are transferred without significant changes to the case of a finite 
closed interval. 

Notice, that ignoring values of discontinuous functions at their points of discontinuity leads to 
a further simplification of the space of distributions. 

2.2. Dynamic functions. In this subsection wc introduce the supplementary notion of a dynamic 
function which is needed for the construction of the space of test function. It generalizes the notion 
of an ordinary function. Let J = [—5,5]- 

Definitionl. A map f : I ^ ^{J) called the dynamic function. 
We denote the value of a dynamic function / at a point t by f{t){-). 

Denote the set of dynamic functions by dF — dF{I). We say that dynamic functions /i, /2 are 
equal, if fi{t){ ) = f2{t){-) in F(J) for all t £ I. Also, the operations of addition, multiplication by 
an element of R and multiplication are defined in d¥ as follows: for dynamic functions /i, /2 G d¥ 
we define their sum /i + /2 G d¥ by the formula 



for all t G /. Thus, dF form an algebra. 

We define the inclusion of the algebra of ordinary functions F into the algebra of dynamic 
functions d¥ by the map: with an ordinary function / G F we associate a dynamic function / G dF 
defined by 



For a given t, we call f(t)(-) a dynamic value of / at point t. A dynamic value of / which is 
identically equal to a constant function is called an ordinary value. If / G dF has an ordinary 
value at a point t, we denote it by f{t) G R. Then a dynamic function is ordinary if and only if 
it has ordinary values at all points. We denote the set of points, where / has ordinary values, by 





/W(-) ^ fit) 



for all tel. 
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U{f)cl. Let D{f)= I \U{f). 

Let / e d¥, f{t) ^ {t t), with dynamic value /(t)(-) = 7(-). Construct a sequence of 
ordinary functions {fn}^=i defined by 

( j{n{t-r)), ifig(^__L,^+_L), 
~ I 0, otherwise 

The sequence {fn}^=i C F is called the sequential representaMon of / G d¥. 

We define the composition o / g dF of a dynamic function / £ dF and an ordinary function 
5 G F by 

(.9°/)W(-)=ff°(/W(-)) 
for all t € I. Accordingly, define the absolute value of / G dF by 

\.f\m-) = \fm-)\ 

for alltGl. Also, we define the support 



supp/ = cl{t e / : /(<)(•) is not ordinary or f{t) ^ 0} 

where cl stands for closure in I. For a given set M c / let us define 

sup / = supsup/(t)(s), inf / = inf inf /(t)(s). 

M teM seJ tei seJ 

We say that / is bounded on M, if supj^^ |/| < oo. 

A dynamic function / is called nonnegative (positive) (denote / ^ (/ > 0)) if f{f){s) ^ 
(/(t)(s) > 0) for all t € I. A dynamic function / is called nonpositive (negative) if — / is 
nonnegative (positive). 

Let f G d¥ he bounded in a right neighborhood of a point t G I. We say that / tends to c G R 
from the right as f — > t+ if for any e > there exists an t] > such that 

(7) sup 1/ - c| < s. 

{T,T+ri) 

We denote the right-sided limit c by f{T+). The left-sided limit and the limit of a dynamic 
functions at a point are defined similarly. Analogously, define the one-sided limits at a, b. 
We say that / G dF is continuous at r G /, if 

(8) /(r) = /(r+)=/(T-). 

Notice that according to the notational convention above the notations used in this definition 
imply that / has an ordinary value /(r) G R. If / is not continuous at r, then we call / discon- 
tinuous at this point. 

Clearly, if in the definitions given above a dynamic function / is an ordinary function, then 
these definitions coincide with the usual ones. 

2.3. Special algebras of dynamic functions. Motivated by our interest in ordinary diff'erential 

equations, we consider the following subalgebras of the algebra of dynamic functions. 

Denote by dG an algebra of dynamic functions / having dynamic values f{t){-) G G(J) and 
such that the one-sided limits 

fia+), f{b-), /(r+), fir-) 

exist for all t G I. We call the elements of dG the dynamic regulated functions (since G is a factor- 
space, dG is also a factor-space; definitions of the supremum and infimum are changed respectively; 
if there is no need, we will not make special remarks, and call the elements of dG the dynamic 
functions). Let us define the norm in the algebra of dynamic regulated functions dG by 

||/IUg = sup|/|. 
I 

We have the following properties of the elements of dG. 
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Lemma 1. Any dynamic regulated function f G dG has ordinary values on I except for at most a 
countable set. 

The proof of Lemmata-lid other statements of Sections 2 and 3 are provided in Section 5. 

Given / e dG, define a function / by /(<) = f{t) for any t e U{f ). We cah / the ordinary part 
of /. As follows from Lemma ^ the ordinary part / is defined everywhere on / except for at most 
a countable set. Denote ord(/) = /. 

Lemma 2. The ordinary part f of a dynamic regulated function f G dG belongs to G, and the 
equalities f{t+) = f{t+), f{t-) = f{t~) hold for all tel. 

Lemma 3. The set of points of discontinuity T{f) of any dynamic regulated function f G dG is 
at most countable. 

Let us define an inclusion of G into dG. To do this we associate with any regulated function 
/ G G a dynamic regulated function / G dG having dynamic values 

(9) -/(r-) for sG [-1/2,0), = /(t+) for s G (0, 1/2] 

for aU tel. 

Let / G dG be such that 

/(t)(-)eAC(J), /(t)(-l/2) = /(t-), /(t)(l/2) = /(t+) 

for all tel. The set of such dynamic regulated functions with operations induced from dG form 
a subalgebra sG C dG. Notice that T{g) = D{g) {g e sG). 

We say that / G sG is a dynamic function of bounded variation, if ord(/) G BV and 

^ varsgj/(i)(s) < oo, 

where D(f) C / is at most countable according to Lemma 12 The set of dynamic functions of 
bounded variation with operations induced from sG forms a subalgebra. The subalgebra of such 
dynamic functions is denoted by sBV and endowed with the norm 

||/IUBV=|/(a+)| + ||/c||BV+ var;(/(t)(.)) 

teD(f) 

where fc G CBV is a continuous part of / = ord(/) G BV (the .Jordan decomposition of functions 
of bounded variation). 

Example 1. The dynamic Heaviside function 6*^ G sBV is defined by 

0^rit) = [ J; III] 0f(r)(.) =/?(■), 

where (3 G AC(J) is such that /5(-l/2) = 0, [3(1/2) = 1. 

E X a m p 1 e 2. Since G C dG, the usual Heaviside function 6*^ G G is in dG. According to the 
definition of the inclusion of G to dG, we have that 

2.4. Test functions. We denote by V the space of classical test functions, i.e., the space of the 
real-valued continuous functions having compact support in / C R, which is endowed with the 
standard topology |Shi84| . 

We denote by T the space of elements (p G dG having compact support supp(^ C /. We say that 
{'/'nlJ^i C T converges to G T if (p„ — > in dG and there exists a closed interval [c, d] C I such 
that supp(/9„ C [c, d] for all n = 1, 2, . . . Thus, the classical space of continuous test functions V is 
contained in T as a subspace. 

Notice that since ior (p e T, g e dG their product g(p belongs to T, algebra T is an ideal in dG. 
Further, if {<^n}$^Li C T and ipn ^ f, then g(pn — > g(p in T. 

We call T the space of dynamic test functions. The proof of the following theorem is similar to 
the proof of an analogous statement for T> |Yos80| . 
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Theorem 1. T is a locally convex topological vector space. 

3. Distributions with dynamic test functions 

Following standard notation |Shi84j . we denote by V the space of distributions with continuous 
test functions, i.e., the space of continuous linear functionals 2? ^ R ,Shi84 . 

We denote by T' the space of distributions with dynamic test functions., i.e., the space of con- 
tinuous linear functionals T — > R. The value of a distribution / G T' on a test function G T is 
denoted by (/, v?) £ R. 

Examples. Let / £ dG (in particular, / e G). We define a distribution, which is also denoted 
by /, by the formula 

(10) {f,^)^ jjm{t)dt 

where £ T ^ (p = ord((/?), / = ord(/). The linearity and continuity of / follows from the 
properties of the integral and the definition of the ordinary part of a dynamic regulated function; 
a distribution / G T' is called a regular distribution. 

Example 4. Define a distribution 5" by the formula 

(11) (5.",^) = jv{r){s)a{s)ds 
where if (^T and a G L( J) is such that 

(12) j a[s)ds = 1. 

The linearity and continuity of (5" follows from the properties of the integral and the definition of 
the convergence in T, so (5" G T'. The distribution 5" G T' is called the delta-function concentrated 
at the point t £ I and having shape a G L(J). Now, for continuous test functions Lp ^T) 

i.^r:V)^ / p{T)a{s)ds^if{T) I a{s)ds ^ ip{t) . 
J J J J 

Thus, the delta-function St G V extends to the family of delta-functions 6" G T' . For a given 
delta-function 6" we construct a sequence 

na{n{t ~ r)), if i G (r - ^, r + ^), 
0, otherwise 

We call {w"}^]^ the delta-sequence having shape a G L(J). 

Examples. Let us define a distribution 6^ E T' , which is also called delta-function, by 

(13) {S^,p) = Xp{T+) + {l-X)p{T-), 

where A G R and the restriction 6^\ti = Sr E T>' . Clearly, (|13|l can not be obtained from (|ll|l . The 
notion of a delta-sequence is not defined for the delta- function (|13|l . 

Remarkl. In the general case, the delta- function in T' is defined as an afhne combination of 
delta-functions (|ll|l and (|13|l . Below we consider mainly the delta- functions of the form (|ll|l . 

In the space of distributions T' the linear operations of addition and multiplication by elements 
of R are introduced in a standard way, so T' is a linear space. Let /« — > / (/«, / G T') in T', 
if ifmP) —> ifi^) as n — > oo for any test function G T. The proof of the following lemma is 
similar to the proof of analogous theorem for V |Shi84l p. 65]. 

Lemma 4. // converges in T' , and ipn ^ ^ in T , then (/„, Lp) — + 0. 
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Theorem 2. Let the sequence {/nj^i C T' he such that for any Lp T the sequence {(/„, <^n)}5^Li 
converges as n —> oo. The functional f defined on T by 

(14) (/,^)- lim(/„,(^) 

n — >oo 

is linear and continuous and thus belongs to T'. 

Theorem 3. Any distribution in V can be extended from V to T. 

Thus, the space T' is a space of extensions of distributions in V from V to T. 

Notice, that the restriction of the distribution ~ 6J (1 T' , where a, 7 are the shapes of delta- 
functions, from T to is a zero distribution in J)' . Hence, since the shapes a, 7 can be chosen 
arbitrarily, addition of 5" — SJ to an extension from P to T of any distribution in V, gives another 
extension. Consequently, any distribution in 2?' has infinitely many extensions from 7) to T. 

A distribution / £ T' is called nonnegative (nonpositive) if for any ip G T such that 1^ > we 
have (/,(/?) > (or {f,(p) < 0, respectively). 

Let us define in T' the operation of multiplication of distributions by the elements of dG. Since 
for any g G dG, ip d T their product gip belongs to T, we may define 

(15) (5/,^)-(/,5^)- 

where f Q T'. In particular, since G C dG, any distribution in T' can be multiplied by any 
piece-wise continuous function. 

Theorem 4. Let fn f in T , gn g in dG. Then gnfn ^ gf in T' . 

The operation of multiplication in T' defined by IjlSfl is continuous, commutative and associative 
in the sense that {gh)f = g{hf) in T' for any g,h dG, f G T'. 

E X a m p 1 e 6. The product of the delta-function (5" and the Heaviside function G sG is given 
by the formula 



(16) {6^,6^, = {5^,ePM = p{r){s)(3{s)a{s)ds. 
If Jj f]{s)a{s)ds 7^ 0, then can be rewritten as 

(17) e^J^^(^j^p{s)a{s)dsyi, 
where 7(-) = /3(-)a(-)/ j j (3{s)a{s)ds satisfies (jEJ. 

Example?. The product of the delta-function 5" and the Heaviside function 6',- G G is given 
by the formula 

(18) (0,(5^?, ^) = {5'^,er^) = ' ^{T){s)a{s)ds. 



1/2 wvn 

If /p a{r)dr ^ 0, then H18|l can be rewritten as 

where the shape of the delta-function 61 is defined by 7(5) — for s G [—1/2,0), 7(3) = 
a{s) I ^l'"^ a{r)dr for s G (0,1/2]. 

R e m a r k 2. Let us show that the equality p7|l can be obtained if (5" and 9^ are replaced by 
the terms of the corresponding delta-sequence {w"}J5^Li and sequential representation {fn}^=i, 
respectively. Let Jj P{s)a{s)ds ^ 0. We have that 

f3{n{t - T))na{n{t - r)), it t e {t - + 
0, otherwise. 



/n(iX(0 

i.e.. 



f3(n{t-Ty)a{n(t-T)) 



fn{t)<it)= / ma{s}ds{ — 

0, otherwise. 



iftG (r- 2^,^+_L), 
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where 7(3) = P{s)a{s)/ Jj a{s)j3{s)ds satisfies ((T^ . 



P{s)a{s)ds ]uj1 



for any n G N. In this sense the operation of multiplication in T' has Property B (see the 
Introduction). 

Examples. In the general case, let us consider the product of the delta-function (5" and a 
dynamic regulated function / G dG: „ 

(/<5^,^) = (<5^,/^)= / ^{T){s)f{T){s)a{s)ds. 



If / has an ordinary value at r, then the result above can be rewritten as 

Remarks. The intersection T n G with the topology induced by T determines a space of 
discontinuous test functions. Since this is a subspace of T, we obtain a factor-space of T', which 
is called the space of distributions with discontinuous test functions (DK05) (see the Introduction; 
also, see IKEMlEirMI)- 

Let g E sWV. Define the derivative g G T' by the formula 

(19) {g,^)^ f md9c{t)+ f ^{r){s){g{T){s)yjs. 

where i^p E gc G CBV is a continuous part of orA{g) G BV. If for any r G T{g) we have 
c^r(ff) = g{T+) - g{T-) ^ 0, then 

(20) 5 = ffc+ 

where 



(21) {9c.v) = jmd9c{t). 

the shapes au are defined by ak{s) — {9{Tk){s))s/<Tr^{g) (s G J, fc G N). 

Theorem 5. The functional defined by U9\} is linear, continuous and thus determines a distribu- 
tion in T'. 

Theorem 6 (Leibniz product rule). For any f,g G sW¥ we have 

(22) if9y =f9 + f9eT'. 
Example9. The derivative of the Heaviside function 9f^ G sBV is given by 

where a = /3 is a shape of the delta-function. 

3.1. Ordinary differential equations with distributions. In what follows, notations 7^' and 
stand for the spaces of n-valued distributions with components in T' and V respectively, 
where convergence, linear operations, operations of multiplication and differentiation are defined 
componentwise. We introduce analogous notations for the spaces of ordinary functions and dynamic 
fimctions BV„, AC„ and sBV„, respectively. 
Consider in 7^ a Cauchy problem 

(23) f{t,x)+g{t,x)S^, x{to~)^xo, 

where G 7, xq G R", and the functions / and g with values in R" and Mn, respectively, are 
continuously differentiable in both variables, (5" G 7^ is a vector-valued delta-function, 

a = (ai, . . . , a„)^, 

the value of (5" on 93 G T is defined by 

(24) {6'^,^)^({S^\^),...,{S^-,^) 
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We define a solution of the Cauchy problem (|23|l to be a dynamic function x G sBV„ which 
satisfies ()23|l in 7^. Define an ordinary solution x ord(a;) G BV„. 

Then there exists a solution of the Cauchy problem (|23|l x G sBV„ having ordinary values x{t) 
for all t ^ T, and having dynamic value x{t){-) (denote 7(-) = x{t){-)), such that 

(25) i(i)=/(i,x(0), 
for alH ^ r, and 

(26) 7(s)=5(T,7(s))a(s), 7(-l/2) = .t(t-). 

(see ExamplelHl, where x(r+) = 7(1/2). 

Notice that in the result of the substitution of a; G sBV„ into H23(l , the operations of differen- 
tiation, composition and multiplication (see example [HJ, which arise in (|23|l . are correctly defined 
in 7^, in contrast to the space V'^ (see the Introduction). 

R e m a r k 4. Consideration of ordinary differential equations with distributions in T)'^ is presented 
in S V97. .MRQS^ .SZQ?, . where, in particular, the following Cauchy problem is considered: 

(27) X ^ f{t,x) + {g{t,x)L)Sr, x{ta)^xa, 

where /, g are the same as above, St G is a scalar delta-function, t = (1, . . . , 1)^. 

As is mentioned in |SZ97| , it is reasonable to expect that in the general case a solution x of (|27|l 
is discontinuous at point t. Thus, in the general case, the notation in (|27|l is incorrect from the 
point of view of distribution theory, since H27|) contains the product of a discontinuous function 
.g(T, x{-)) and a distribution dr G V, which is undefined in V'^. In |SV97[IMn96llSZ97j the following 
definition of a solution of (|27|l is proposed. The left-continuous function x G BV^ is said to be a 
solution of lf?7|) . if there exists a delta-sequence {ujm}m=i: — > in 2^', such that 

(28) Xjn X 

(m —>■ go) in a weak topology of BV„, where Xm G AC„ is a solution of the approximating problem 



(29) X = f{t, x) + {g{t, x)L)ujrn{t), x{to) = a;o. 

The necessary and sufficient condition for the uniqueness of solution of H27|) . i.e., its independence 
on the choice of {uJm}'i^^i, is a Frobenius condition 

i—1 2—1 

for aU t G /, X G R", fc, j, m = 1, . . . , n, 5 = (5y)y-=i |SZ97[ lMii93| . If ^ is satisfied, then there 
exists a unique solution x G BV„ of (I27II . such that 

(31) i = /(t,x(0) 
for alH ^ T, and 

(32) Ti{s) ^ g{T,ij{s))i, r,(-l/2)-x(T-), 

where 2;(r-|-) = r/(l/2) [^tT???! . 

Notice that (|32|l coincides with H2t)|l for a = i. Indeed, for H2t)|l the Frobenius condition (|30|l is 
a necessary and sufficient condition for the independence of the value 7(1/2) on the choice of a 
satisfying (|12|l |Cuy83| . Also, if (|30|l is satisfied, then a solution of (|27|l coincides with the ordinary 
solution of ((221) • 

The comparison of the results for equation (|23|l , which is considered in the space of distributions 
7^, and the results for equation H27|) . which is considered in the classical space of distributions 2?^, 
shows, that in contrast to the approach based on the space V'^ |S V97I IMR96I [S797] . the use of the 
space of distributions 7^ allows us to make notations correct from the point of view of distribution 
theory, and to eliminate the restrictive Frobenius condition l|3U|) . which can not be satisfied in 
many problems of optimal control theory |Mil93| . 
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4. Multiplication of distributions by discontinuous functions in other approaches 



4.1. Multiplication by discontinuous functions in the Schwartz space of distributions. 

Various definitions of tlie product of a distribution and a discontinuous function in T)' were proposed 
in SZ97. .Fil88 KBm <Kur96[ ITvr02| . and also in ,Sar95l ISar94[ ISar03| . wliere tlie operation 
of multiplication of a distribution by a discontinuous function is induced by the operation of 
multiplication of two distributions. 

Furthermore, we concentrate on the multiplication of the Heaviside function 9^ and the delta- 
function. The product OrSr e V is defined in |SZ97l iFiISSl IKB98I IKutMI ITvrn2l ISar9Rl ISar94[ 
ISar03| by the formula 

(33) OrSr = c5r, 

where c G R is different in different definitions. In particular, in |Sar95l rSar94l [Sar03j . where 
the family of the products of distributions {■a}ae'D is defined (here I? is a space of infinitely 
differentiable test functions), for any G P we have 

(34) Or 6r = <7{0)Sr. 

The value c = 1/2 is considered in |Fil88l IKB98I IKur96j V If c = 1/2, then we have that 
[O-rOr)' — OrSr + S^Or, i.e., Lcibuiz product rule holds. Also, the choice of the value c = 1/2 is 
motivated by the symmetry of the delta- function 6r € V. Another definition of the form (|33() was 
proposed in |SZ97j , so the solution of the Cauchy problem H27|l in coincides with the solution 
(128(1 . In |SZ97j the value of c G R is determined by the right-hand side of the equation 127|) . 

Nevertheless, any definition of the product of the form H33() leads to contradictions with distri- 
bution theory, which are unavoidable in V . Namely, let c G R be given. Then we have to be able 
to replace delta-function in H33|l by the terms of arbitrary delta-sequences {oJn}^=i, in V 

so that 

(35) lim OrLOn ~ c lim A„ 

n — ^oo 71 — ^oo 

in V . Nevertheless, for a given c G R one can easily find delta-sequences {wn}5^Li> {-^n}$^Li such 
that ((35|l is not true. This implies that the operation of multiplication ((33|l is not continuous in 
V . Also, if c / 0, c 7^ 1, then the operation of multiplication 1)33(1 is not associative. 

In T', the product of the Heaviside function 9r and the delta-function 6" G T' is given by the 
formula 

(36) 9r6^ = (^J^ ' a{s)ds^6:i, 

(Example Ej), where S'J G T'. Let us notice the following: 

1) Any definition of the form 1(33(1 can be obtained from ((36(1 for the delta- functions 6" G T' 
having the shape a such that , 

a{s)ds — c. 





2) In contrast to ((33(1 . the operation of multiplication ((36() is associative, satisfies the Leibniz 
product rule (ThcorcmEJ, continuous and, furthermore, has Property B (see the Introduction and 
Remark EJ. 

Thus, consideration of the product of the Heaviside function and the delta-function in T' allows 
us to avoid the aforementioned contradictions with distribution theory. 

4.2. Multiplication by discontinuous functions in the Colombeau generalized functions 
algebra. As is mentioned in the Introduction, the general multiplication of distributions is impos- 
sible in the classical space of Schwartz distributions with continuous test functions. This defect of 
the Schwartz distribution theory lead to the construction of the Colombeau generalized functions 
algebra G, containing the space of distributions V as a subspace (here V is a space of distributions 
with infinitely differentiable test functions), where such general multiplication exists ^Col85| (see 
description of the Colombeau algebra in |Col92l IE.()97p . The algebra of Colombeau generahzed 
functions was applied successfully to many problems related to ordinary and partial differential 
equations, e.g., sec ,H097..CLNP89...CH92..,Col92, . In the classical approach the Heaviside function 
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Ot (z G and the Dirac delta- function 5^ ^ Q are defined as the derivatives of continuous functions 
in so Or, 5t ^T)' (Z Q and the product 

eg\v' 

is not a distribution |Col85j . However, as was observed in jCLNP89j .^ instead of consideration of 
the elements of the canonical embedding V C Q, often it is natural to work with the algebra Q 
itself. Namely, a generalized function f £ G is said to have 9t or St as its macroscopic profile, if 

f KiOr 01- f ~ Sr, 

respectively, where w stands for the association in Q. If / has 0t- or 6t as a macroscopic profile, 
then / is called the Heaviside generalized function or the Dirac generalized function, respectively 
|CLNP89] . According to |CLNP89| , there exist different Heaviside generalized functions and Dirac 
generalized functions, which differ by their microscopic profiles. The equalities of the form 

(37) Oi^rS2,T = cSs^r, 

where 0i^r £ Q is a Heaviside generalized function, S2^t, "^s.r G G are two Dirac generalized functions, 
were studied in |Col92| . where the value of constant c e R was derived from the supplementary 
partial differential equation, arising from the statement of the physical problem. As follows from the 
definition of the algebra G |H097| . equality (|37|l can be also obtained if the Heaviside generalized 
function and the Dirac generalized functions in l|37(l are approximated by continuous functions. 

The equality H37|) is similar to the equality Hl()|l in the space of distributions T', where the 
constant c G R is given explicitly (see also Remark |5J). However, despite this similarity, space T' 
also possesses other properties. In particular, T' contains delta- functions H13() (see Remark QJ, 
which allow us to unify in the subsequent paper different definitions of solution of the ordinary 
differential equations with distributions (see |Fil88| on ambiguities arising in the definition of the 
notion of solution of the differential equations with distributions). 



5. Proofs of statements of Sections 2 and 3 

Proof of Lemma^ Suppose that there exists an uncountable set A <Z I such that f{t){-) is not 
ordinary for all t £ A. Let us define a function g : / ^ R by 

(38) 5(i)=sup/(t)(s)-inf /(t)(s) 

seJ 

Since / G dG is bounded, ^ g{t) < oo for all t £ I. Obviously, / has an ordinary value at t if 
and only if g{t) = 0, so g{t) > for all t £ A. We define 

A„^\teA: ^— < g{t) ^ -\ {n ^ 1,2, . . .), Ao ^ {t £ A : git) > 1}. 

Then A = U^^gAn- Since A is uncountable, there exists no ^ such that Ang is uncountable. 
Without loss of generality we may assume that no ^ 0. Since An^ is uncountable, there exists 
a closed interval [c, d] C / such that A'^^^^ = Ang n [c, d] is uncountable. Consequently, due to 
compactness of [c, d] , the set A'^^ has a limit point ^ £ [c, d] . Without loss of generality we may 
assume that there exists a right-sided neighborhood of ^ which contains infinitely many points 
of A'^^. Let us show that this implies that f{£,+) does not exist. Suppose that p — f{£,+) £ R 
exists. Then by definition for any e > there exists 77 > such that ||/ — < e, i.e., 

suptg(a,6sup^gj \f{t){s)-p\ < e. Now consider the value of sup^gj 1/(0 («) for t e (^,C + '7)n 
A'ng , where (^, i + r])n A'^g ^ 0. Then the inequality 

(39) sup|/(t)(,s) -p\ ^ Vsup(/(i)(s)) - inf (/(t)(s))) 



We would like to thank Professor J.F.Colombeau who kindly mentioned these results to us. 
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holds (see proof below). As follows from the definitions of ^'^^ and g, 

i(sup(/(i)(s))-inf(/(t)(s)))^l/no 

for aU t e A'^ig n {(,,£, + 77). Consequently, suptg(^ sup^^j \ f{t){s) - p\ > Since for any 

7^ > A'jjj n , ^ + 77) 7^ 0, the last inequality holds for any 77 > 0. This contradicts to the definition 
of the right-sided hmit p — /{£,+)■ 

Let us show that the inequality holds. If f{t){s) —p ^ for all s G J, then sup^gj \f{t){s) — 
p\ — sup5gj(/(<)(s) — p) and holds, where J^^j{f{t){s) — p) > 0. In the general case we can 
represent f(t)(s) — (/^(s) — 9 (s), where q'^^s), q^{s) ^ 0. Then 

sup \f{t){s) ~ p\ — max{sup(7^(s), supg^(s)} ^ l/2(sup(7+(s) + sup(7^(s)). 

sej s£j sej seJ sej 

Since sup^g j(/(t)(s) - p) = sup^^jq+{s), - mise.j{f{t){s) - p) = sup^gjg~(s), we obtain that 
(EH is true. □ 

Proof of Lemma\^ For a given r G / let us denote p = /(r+). Then for any e > there exits rj > 
such that sup(^ ^._|_,^) \p — f\ < £• According to Lemnia^ the intersection (r, r + 77) n [/(/) 7^ 0. 
Then 

sup \p-f{t)\= sup sup|p-/(t)|^ sup \p-f\ 
te{T,T+ri)nus{f) te{T,T+ri)nus{f) seJ (T,r+?))nus(/) 

where the first equality holds since \p — f{t) \ does not depend on s G J. Thus, for any t G (t,t + 77), 
where f{t) is defined, i.e, for all t G (r, r + 77) nU{f), we have |p — f{t)\ < s. Then the right-sided 
limit f{T+) exists, which is equal to /(t+). The proof for the left-sided limit is analogous. Since 
T G / was chosen arbitrarily, the statements of the lemma holds. □ 

Proof of Lemma\^ The proof follows from an analogous statement for the space of regulated func- 
tions G |Hon75j . the definition of a point of discontinuity of a dynamic function, and Lemma 

El □ 

Proof of Theorem\^ The linearity of the limit functional is obvious. Without loss of generality it 
suffices to show that / is continuous at (p — only. Let (p„ in T. Suppose that (/, tp „)-/->■ 
(n 00). Being considering, if necessary, a subsequence of {{f,^n)}^=i, we may assume that 
for any 77 = 1,2,... the inequality |(/, > Eq holds for certain Eq > 0. Due to p4(l for 
any k — 1,2,... there exists 77^,. such that Hfn^jfk)] > Without loss of generality we may 
assume that = k, i.e., |(/fc, fk)\ > ^ for any fc = 1, 2, . . . The last inequality contradicts to the 
conditions of the lemma. Thus, (/, ipn) —>■ asn 00, i.e., the limit functional / is continuous. □ 

Proof of Theorem\^ According to Theorem ^ the space of dynamic test functions T is locally- 
convex; I? is a subspace of T. Then according to the Hahn-Banach theorem |KA82| , every linear 
continuous functional defined on a subspace V has an extension to the whole space T. □ 

Proof of Theorem^ Note that g„</j — > gip in T for any ip ^ T. Consequently 

\{9nfn,ip) - (g/, V)l = \{fn,gnip) " < " (/«,5</')l + 

+ \{fn,g<f) - {f,g<f)\ ^ \{fn,gnf - gf)] + \{fn,g(p) - {f,g'p)\ 

(the first and the second summands tend to zero according to Lemma 0] and due to convergence 
/„ ^ / in T', respectively). □ 

Proof of Theorem^ The Stieltjes integral in the right-hand side of p9|l exists since (p ^ gc ^ 
(see jDern2| V and the convergence of the series follows from the inequality 

f{T)[s){9{T){s)yjs sup|(y9(T)(s)| / |5(t) (s);, | ds sup | (^(t) (s) | var.g J (^(t) (s)) 

sGJ J J sGJ 

for all r G T{g). Then the following inequality holds 

Y Jj^i'^)(^)i9{T){s))sds ^sup|<y9| Y var,g./(5(T)(s)), 



reT{g} 
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where convergence of series in the right-hand side follows from the definition of the algebra sWV. 
Thus, the value of g is defined for all test functions (p ^ T. The linearity and continuity of g 
follow from definition of the convergence in T, properties of the Stieltjes integral and properties of 
convergent series. □ 

Proof of Theorem\^ We have 

(40) {fg,^) = if,g^)= f mmdfcit)+ ^ /'g(r)(sMr)(s)(/(r)(s)),ds, 

(41) {fg,^)^igj^)= f mfit)dgc{t)+ J2 //(T)(s)^(r)(s)(5(r)(s)),ds. 

In l|4U|l and H41|l we can perform the summation by T {f)UT (g) , since we need to add zero summands 
only. Consequently, 

(42) [ 9ir){sMT){s)ifiT){s)y,ds+ / /(r)(s)^(r)(s)(.9(r)(.)),rfs = 

= E U{r){s){g{r){s)f{T){s)y,ds= I ^{r){s){g{T){s)f{r){s)y,ds, 

where the last equality is due to T{fg) C T(/) U T{g) and r e T(/) U T{g) \ T{fg) if and 
only if /(r) = or (/(r) = 0, so we may exclude in (|42|l the summands corresponding to t G 
T{f)UTig)\T{fg). We have 



(43) / m9mc{t)+ / mf{t)dAt)= / md[ / m)dfc{o+ / fiodgdo], 

Jl Jl Jl \Ja Ja 

where I — {a,b). Let us show that the following equality holds 

(44) ifgUt) = f giOdfciO + f KOdgdO + f{a+)g{a+) 



for all t € I. Now {fg)c = {f9)c, fc = fc, 9c ~ gc, thus we may prove H44I) for the ordinary parts 
only. Hence 

f{t)g{t) = (/c(t) + /;(t))(5c(t) +5/.(i)) = /c(i)ffc(i) +/c(t)5/.(i) +3c(t)A(i) + A(t).9/.(t) 
for all t e /. Thus, 



tmdm)+ f fiOdgc 

J a J a 



(0 + /(a+).g(a+) 



9cmfc{i) + / 9h{mc{i) + / fcmui) + / jhmm) = 

J a J a J a 

= fc{t)9c{i) - fc{a)9c{a) + f{a+)g{a+) + fc{t)gh{t)~ 

- E fM){9{n+)-g{n^))+mh{t)^ E 9c{n){f{n+)-f{n-)), 

a<ri<it a<Ti<t 

for all i G /, where fc{o)fc{o) = f{a-\-)g{a-\-) and 

fc{t)gh{t)- E fc{n){9{n+) - 9{n-)) 

a<Ti<t 

is a continuous part of /cffh, 

9cit)fh{t)- E .9c(t.)(/(t,+) - /(r,-)) 

a<Ti<t 

is a continuous part of gc//i, so equality (|44|l is true. The comparison of (|4()|l - (|44|l with the 
definition of the derivative in T' gives us (|22|- CD 
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